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ABSTRACT 

Models are presented for the Lya emission of dusty high-redshift galaxies by 
combining the Press-Schechter formalism with a treatment of the inhomogeneous 
dust distribution inside galaxies. It is found that the amount of Lya radiation 
escaping from the galaxies strongly depends on the time over which the dust 
is produced through stellar activity, and on the ambient inhomogeneity of the 
HII regions that surround the ionizing OB stars. Good agreement is found with 
recent observations, as well as previous non-detections. Our models indicate 
that the dust content builds up in no more than ~5x 10 8 yr, the galactic HII 
regions are inhomogeneous with a cloud covering factor of order unity, and the 
overall star formation efficiency is at least ~ 5%. It is predicted that future 
observations can detect these Lya galaxies upto redshifts of ~ 8. 

Subject headings: cosmology:theory - early universe - galaxies:evolution - 
galaxies:star formation - galaxies:ISM 
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1. Introduction 

One expects that the ionizing radiation from young stars leads to prominent Lya 
emission through the recombination of hydrogen in the interstellar medium (Meier 1976). 
Because the Lya line is narrow and strong, it stands out against the continuum background 
and should provide a signature of primeval high-redshift galaxies (Partridge & Peebles 
1967). Until a decade ago, the search for these high-redshift Lya galaxies had enjoyed no 
compelling successes (e.g. Djorgovski 1992; Thompson et al. 1995). A large class of high 
redshift radio galaxies was discovered in recent years (of. Chambers & Miley 1990), but 
their numbers are still significantly below the theoretical predictions of about 10 3 ~ 5 of such 
primeval galaxies per square degree (cf. Pritchet 1994). As was pointed out in the above 
paper, the lack of detections was becoming a source of confusion and concern. However, 
in the last couple of years, with improved sensitivity on large area telescopes, these young 
galaxies are finally being detected (e.g. Hu, Cowie, & McMahon 1998; and references 
therein). Clearly, this population is of great interest to the field of galaxy formation and 
the early evolution of the universe. It is mandatory to understand why earlier surveys have 
been unsuccessful, how many objects are still expected to be found, and what physical 
conditions pertain in these early systems so that the Lya radiation may escape. 

On galactic scales, the spatial distribution and the amount of stellar dust are crucial 
ingredients for the radiative transfer effects in the Lyct line, such as resonant scattering 
and dust attenuation, and will determine the emerging Lya luminosity (e.g. Neufeld 1991). 
Indeed, for a homogeneous dusty HII region it is well known that only a negligible fraction 
of the Lya photons can escape the medium (of. Spitzer 1978). Observations by e.g. 
Hartmann, Huchra, & Geller (1984) and Terlevich et al. (1993) have firmly established the 
strong decrease in Lja equivalent width with increasing oxygen abundance. 

The presence of dust in galactic HII regions therefore needs to be studied in detail. 
In a cosmological context, the formation of dust requires the presence of metals, and is 
therefore intimately related to the star formation process itself (Miralda-Escude & Rees 
1998). Connected to this metal production is the question of the early enrichment of the 
intergalactic medium (IGM) as seen in damped Lya systems, as well as that of the hot 
cluster gas observed in lower redshift systems. The interplay between the IGM and the 
galactic environment likely plays an important role in the distribution of metals and dust 
through galactic wind expulsion and accretion (Ferrara 1998). These effects are particularly 
important in the light of recent observations by Pettini et al. (1998). Although their 
signal-to-noise ratio should be improved, these observations suggest that damped Lja 
systems do not trace compact dense regions, i.e. primordial galaxies, but reflect the more 
diffuse IGM enriched by these same primeval structures. As such, the metal content of a 
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galaxy or the IGM may show two quite different sides of the same star formation process. 

The importance of dust has been encountered in the past both observationally and 
theoretically. Obscuration effects have been studied for quasars by Fall & Pei (1993) and 
galaxies by Chariot & Fall (1993). The latter authors perform a study somewhat similar 
to the present paper, but here we have included the effects of inhomogeneity and large 
amounts of dust, and also utilized the physically motivated Press-Schechter formalism for 
the star-formation history. Local starbursts are known to emit the bulk of their energy 
in the infrared through dust absorption and re-emission. The evolution of the universal 
average star-formation rate (Madau et al. 1996) must likely be modified for dust extinction 
to reflect the true star formation history of the universe. Finally, recent SCUBA (JCMT) 
results indicate the presence of extremely dusty galaxies (Smail et al. 1998) at high redshifts 
(z ~ 3), which could also be related to the production of the recently discovered cosmic 
infrared background radiation (Puget et al. 1996; Hauser et al. 1998). 

The aim of this paper is to investigate the importance of dust for the recently discovered 
Lya emitters; to show that it is possible to derive constraints from the observational data 
on the formation of galactic structure, and on the first generation of stars; and finally 
to make predictions for higher, previously unexplored redshifts. It is assumed that the 
bulk of the Lya photon production is driven by stars, and that the contribution of type 
II supernova shocks or non-thermal emission is negligible (Chariot & Fall 1993). An issue 
which will not be addressed here in detail is the formation of Lya galaxies (e.g. Primack 
et al. 1998). Although this is clearly an important question, the aim here is to set up a 
general cosmological model, based on the Press-Schechter formalism, where the absolute 
star formation rate and the physical state of the system are varied while the observational 
consequences are explored. 

2. Lya Emission from HII Regions 

It is well known that even a modest amount of dust inside an HII region is sufficient 
to attenuate all of the produced Lya radiation, because the line optical depth is of the 
order of 10 4 . However, this situation can be alleviated if the dust content is negligible, or, 
more interestingly, if the medium is inhomogeneous (Neufeld 1991). The latter author has 
argued that in a multi-phase medium the escape of Lya radiation is significantly enhanced 
compared to the standard case of resonant scattering in a homogeneous medium. 

The radiative transfer problem which needs to be solved for the Lya line is well 
studied and has been investigated extensively (e.g. Adams 1972; Hummer & Kunasz 1980; 
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Neufeld 1990). Our aim here is to model an individual galaxy with a range of masses for 
the ionizing stars, dust content, and inhomogeneity. A numerical Monte Carlo approach 
is therefore adopted whose general implementation can be found in Spaans (1996; and 
references therein). In these computations, a Scalo IMF is assumed for the spectral types of 
the central stars in the HII regions, ranging from 05 to Bl. The stars are distributed in a 
statistically homogeneous manner inside a percolating multi-phase medium. A typical line 
width AV of 8 km/s has been used. Care was taken to sample the line profile sufficiently 
far into the wings, more than 10 2 (A r 2oAVs) 1//3 km s _1 for a homogeneous medium with a 
hydrogen column density N 2 o in units of 10 20 cm" 2 and a velocity dispersion AV5 in units 
of 10 5 cm s _1 (Harrington 1973). For an inhomogeneous medium the formalism of Neufeld 
(1991) for a multi-phase medium is adopted. 

Analogously to the Neufeld investigation, we parameterize the multi-phase medium 
by opaque clumps embedded in an inter-clump medium of negligible opacity. The clump 
covering factor F cov , i.e. the average number of clumps along a line of sight, then fixes the 
degree of inhomogeneity. For the work presented here, the escape fractions F esc = exp(— r#) 
were computed on a grid of models with various dust contents, Z d = 1CT 2 — 1 solar (Draine 
& Lee 1984), and covering factors, F cov = 1, 5, 00. The results of the numerical computations 
for these three covering factors are shown in Figure [1|. As this figure shows, the covering 
factor has a significant effect on the escape fraction. In particular, the inhomogeneous 
percolating slabs are much more transparent than the homogeneous ones; the difference 
around Z d ~ 10" 1 solar is over an order of magnitude. The reason for the large difference is 
that photons incident upon highly opaque clumps will be reflected back into the interclump 
medium. The photons do not penetrate the opaque clumps and spend most of their time in 
the interclump medium, where the opacity is very small, facilitating their escape through 
repeated reflections off the clumps (Neufeld 1991). 

In order to check the accuracy of our computations, we performed similar calculations 
for the slab geometry adopted by Neufeld (1991). We were able to obtain good agreement 
with the analytical solutions for these geometries. In particular, the result for the effective 
optical depth r* « (AV/AV*) 1/2 with 

AV. = 680F COT 2 ^ 1/2 x(H°), (1) 

was reproduced with high accuracy for the neutral hydrogen abundance x. For 
computational convenience, the numerical calculations were terminated once the escape 
fraction reached 10~ 2 , although the results presented below do not depend on this choice. In 
all of the calculations below, we assume the medium to be homogeneous until its metallicity 
exceeds 3% solar, and allow for inhomogeneities once this level of enrichment is reached. 
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3. Cosmological Abundance of Lya Emitters 

In order to model the cosmological abundance of high-redshift Lja emitters as 
a function of redshift, we assume that the formation of dark matter halos follows the 
Press-Schechter (1974) theory. Hence, the net rate of change in the comoving number 
density of halos with mass Mh a i is given by the derivative ^(dn c /dM), where dn c /dM is 
the Press-Schechter mass function (comoving number density per unit mass). This rate 
includes a negative contribution from the disappearance of small halos in merger events; 
the net -^{dn c / dM) becomes negative when mergers dominate, i.e. for masses below a 
characteristic mass scale M*(z) at any given redshift z. We make the simplifying assumption 
that no halos form with masses below Mh a i < M*(z), and therefore set -j^(dn c / dM) = for 
halos below this mass. 

Next, we assume that every halo forms a galaxy that goes through a Lja emitting 
phase. Indeed, after a halo collapses and virializes, its gas can fragment into stars, 
provided it can cool efficiently below its initial virial temperature. Because of the lack 
of a significant amount of H2 molecules at high redshifts (Haiman, Rees & Loeb 1997), 
the leading cooling mechanism is expected to be collisional excitation of atomic hydrogen. 
The corresponding requirement for halos to cool efficiently is that their mass is at least 
M min ~ 1O 8 M [(1 + z)/ll]- 3/2 (Haiman & Loeb 1998). We therefore assume that in halos 
that exceed this mass, a fraction e* of the gas turns into stars, with a constant rate of 
star formation over a period of t* years. The result of this process is a dwarf galaxy 
with a stellar mass of M star = e*(f2 b /fi )M ha i o . The values of and could depend on 
several further parameters, such as the formation redshift, halo mass, or the initial angular 
momentum of the gas. There could be further complications, such as a time-dependent star 
formation rate. A treatment of these important issues is beyond the scope of this paper, 
and for simplicity we adopt the simplest assumption, i.e. that e* and t* both have the same 
constant values in each halo. 

Finally, an important effect we must include is an external feedback from photo- 
ionization by the UV background, which heats the gas before it is able to cool and condense 
inside the dark matter potential wells. This feedback is an inevitable consequence of the 
background UV flux that builds up after the re-ionization epoch, when the cosmological 
HII regions have overlapped and the universe is transparent to the ionizing flux from each 
individual source. Several authors have discussed the consequences of this feedback from 
first principles (Babul & Rees 1992; Efstathiou 1992; Quinn, Katz & Efstathiou 1996; 
Thoul & Weinberg 1996; Navarro & Steinmetz 1997), and concluded that the collapse of 
gas is inhibited in halos with circular velocities below t> circ = 30 — 80 km s -1 . The apparent 
lack of high-redshift quasar candidates in the Hubble Deep Field (HDF) also requires the 
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existence of some type of feedback of a similar magnitude. To be consistent with the HDF 
data, Haiman, Madau & Loeb (1998) have found that either the formation, or the fueling 
of black holes was prevented in halos with t> circ < 75 km s" 1 . Accordingly, we impose a 
minimum circular velocity t> C i rc = 75 km s _1 for the halos of luminous galaxies that form 
after the re-ionization epoch, here taken to be z re i on = 10 (cf. Haiman & Loeb 1997). This 
constraint results in a minimum halo mass of M min ~ 1O 1O M [(1 + ^)/ll]- 3 / 2 , two orders of 
magnitude larger than the minimum mass obtained from the cooling argument above. 

We note that recently discovered high-redshift galaxies (Steidel et al. 1998) have 
revealed Lya emission only in a fraction of all the sources. This may be consistent with 
the small duty-cycle of the Lya emitting phase considered here; alternatively, could 
be increased somewhat if only a fraction / of all halos would undergo a Lja emitting 
phase. In order to keep the predicted number density constant, t* would then need to 
be increased approximately to ~ £*//■ The Press-Schechter halo formation rate also 
depends on the cosmology and power spectrum. For our "standard model" we have 
adopted the concordance model of Ostriker & Steinhardt (1995), i.e. a flat ACDM model 
with a slightly tilted power spectrum (Q , Q^, f^, h, cr^-i, n) = (0.35, 0.65, 0.04, 0.65, 0.87, 
0.96). Convenient expressions for the differential volume element, luminosity distance, and 
time-redshift relation in this model are given in terms of elliptic integrals by Eisenstein 
(1997); a useful fitting formula for the growth function is given by Carroll et al. (1992). In 
addition to this A cosmology, below we consider an open model with the same parameters, 
except the power spectrum normalization is then changed to a sh -i — 1-15. 



4. Model Description 

In order to compare our models directly with observations, we need to compute the 
surface density of Lya emitters on the sky above a given observational threshold flux. As 
discussed above, each of our model galaxies undergoes continuous star formation for 
years. During this time-interval, we relate the local star formation rate (SFR) to an intrinsic 
Lya luminosity, assuming case B recombination, and using Kennicutt's (1983) relation 
between the SFR and Lya luminosity, resulting in L a ^ nt = (M/M yr -1 ) 10 42 erg s -1 . Here 
M = fib/^o-^haio e *A* is the adopted star formation rate. We assume that each model 
galaxy produces stellar dust during the same time-interval, and that some fraction x g of 
this dust stays within the galaxy and enriches its interstellar medium, while the remaining 
fraction 1 — x g is blown out of the galaxy by supernovae, stellar winds, or as a result of 
mergers between galaxies (Gnedin 1998). Galaxies that form subsequently will be assembled 
from an intergalactic gas that is pre-enriched with this ejected dust. 
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It is difficult to estimate from first principles either the fraction x g , or the size 
of the region around each galaxy polluted by the blown-out dust (although see recent 
3-D numerical simulations by MacLow & Ferrara 1998). Instead, here we adopt a 
phenomenological approach, and introduce two additional parameters that characterize the 
overall dust content of galaxies. The first parameter, Z^isuit) is the mass fraction of dust 
(in solar units) in the interstellar medium of each galaxy, due only to the self-enrichment 
from the retained dust fraction x g . Consistent with our assumption of a constant star 
formation rate, we assume that the amount of retained dust increases linearly with time for 
years, until it reaches the final value of -ZdjsM^*)- Similarly, we assume that each galaxy 
deposits dust into the surrounding IGM at a constant rate for t± years, and introduce 
a second parameter, Z d ^ GM (z), denoting the resulting redshift-dependent average dust 
content of the polluted regions within the intergalactic medium. Note that the ejected dust 
may not be fully mixed into the IGM, and the universal average dust fraction of the IGM 
could be lower than Z^igm- Here we avoid this issue by defining Zd.iGM^) to denote the 
average dust content only within the pockets of intergalactic gas from which subsequent 
galaxies form. Observationally, the metallicity of gas in galaxy clusters upto redshift z = 1 
is approximately 1/3 solar, which we take as the representative value for Z d j GM (z = 1). 
Figure || shows the evolution of Z d) i G u(z) in our models, obtained by summing over the 
dust output of all Press-Schechter halos, and assuming three different values of i*. Note 
that the longer the dust-producing phase in each individual galaxy, the steeper the buildup 
to the average Z d j GM (z = 1) has to be, resulting in a lower dust content at higher redshifts. 

In summary, a galaxy that forms at redshift Zf has, at a later redshift z, a total dust 
content of Z^iqm(z) + Z dt isM{At), where At is the time elapsed from zj to z. Both types of 
dust contribute to reduce the intrinsic Lya luminosity L a ^ int implied by the star-formation 
rate. This reduction is expressed by the escape fraction as discussed in Section 2 above. 
Apart from the underlying cosmology, our model therefore has five adjustable parameters, 
F cov , e*, ZdjGM^ = 1), and Zd.istvi^*), that uniquely determine the number density of 
Lya emitters at any flux and redshift. More precisely, at redshift z, the comoving number 
density n c (z, F a ) of emitters with observed line flux between F a and F a + dF a is given by a 
sum over halos of different formation redshifts and ages that exist at redshift z, 



where the factor dM^io/dFa converts the number density per unit mass interval 
to number density per unit line flux interval. The line intensity F a is given by 
F a = F esc (z, z')L a ^ nt (M hll i , e*, t*) / A.irdi,(z) 2 , where d^z) is the standard cosmological 
luminosity distance. 
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5. Results and Discussion 

We first define a "standard model" by the set of parameter values t+ = 5 x 10 8 yrs, 
Fcov — 5, e* = 10%, and Z d j GM (z = 1) = Zci,ism(£*) — 0.3 Z . The first three of these 
values are broadly consistent with the star formation rates, covering factors, and stellar 
mass fractions estimated in low-redshift galaxies; the latter two are chosen based on the 
metallicity of gas in galaxy clusters and within the interstellar medium of low-redshift 
galaxies (c.f. Mushotzky & Loewenstein 1997; Lada, Evans, & Falgarone 1997; Ho, 
Filippenko, & Sargent 1997; Hammer et al. 1997; Young et al. 1996). In Table [l] we 
summarize the parameters of our standard model, as well the ranges we have considered for 
each parameter. 

In Figure ||| we show the resulting surface density of Lya emitters in our standard 
model, with fluxes above different values of the detection threshold. We have chosen 
the quoted 5a detection limit of the recent narrow-band survey of Hu et al. (1998), 
F = 1.5 x 10~ 17 erg cm' 1 s _1 , as a fiducial threshold value. We also indicate the results 
from this survey for the surface density of emitters at the two redshifts z = 3.4 and z = 4.5. 
As the figure shows, the model is in good agreement with these two available data points. 
Although the agreement is encouraging, it is important to examine the sensitivity of this 
result to changes in the model parameters. In Figure [3], we demonstrate the effect of 
changing the star formation rate by changing t*, while leaving e* fixed. This corresponds 
to changing the SFR but not the total stellar mass. The dashed lines show that when t* 
is decreased by a factor of 10, then both starlight and dust is produced earlier, i.e. sooner 
after the formation of each galaxy. This causes the abundance of emitters to peak at higher 
redshifts, but this makes the agreement with the data points only slightly worse. On the 
other hand, when t+ is increased by a factor of 10, the abundance of emitters peaks at lower 
redshifts, and the discrepancy with the data points becomes significant: the model now 

Table 1: The assumed parameter values in our standard model and its variants. 



Parameter 


Standard Model 


Range Considered 


U 


5 x 10 8 yrs 


5 x 10 7 ~ 9 yrs 


± cov 


5 


1 — oo 


e* 


10% 


2 - 20% 


Zd,lGM.{z = 1) 


0.3 Z 


0.1 - 1 Z 




0.3 Z 


0.1 - 1 Z 
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underpredicts the abundance by more then an order of magnitude at z — 4.5. Finally, we 
show in Figure ^| the effect of changing the clump covering factor -F cov . The dotted lines 
show that when the medium is assumed to be inhomogeneous with F cov = 1, the model still 
fits the abundance at z — 4.5, but overpredicts it by a factor of ~ 2 at z = 3.4. On the other 
hand, for a homogeneous medium (F cov = oo), the escape fraction drops dramatically (cf. 
Fig. []]), resulting in a substantial decrease in the surface density of emitters, again leading 
to a significant (two orders of magnitude) discrepancy with the Hu et al. data at z — 3.4. 

It is interesting to note that the surface density of emitters is a steep function of the 
detection threshold. In particular, Figure [3] reveals that, depending on the values of the 
other parameters, raising the flux threshold by an order of magnitude can decrease the 
predicted number of emitters by 2-3 orders of magnitude in the relevant redshift range, 
z = 2 — 5. This is an important feature of the present models, when we consider the existing 
data together with the lack of detections of Lya emitters in previous surveys. In particular, 
based on a lack of detections, Thompson, Djorgovski & Trauger (1995) have deduced a 
limit on the comoving volume density of emitters that is a steep function of the detection 
threshold. The surface density found by Hu et al. (1998) is barely consistent with this limit 
at the faint end of the threshold (LL yQ ~ 10 42 erg s _1 ). Consistency with the Thompson, 
Djorgovski & Trauger (1995) limits then requires that the model predicts a decline in the 
surface density with increasing flux threshold, and that this decline is at least as steep as 
the inverse square of the threshold. As Figure |3] shows, our models indeed have this sharply 
declining feature in the redshift range 2.78 < z < 4.89 of Thompson, Djorgovski & Trauger 
(1995), and can therefore simultaneously explain the new detections without violating the 
existing upper limits on the abundance of emitters. 

Finally, in Figure ||, we demonstrate the effects of changing the cosmology, the star 
formation efficiency e*, or the intergalactic and galactic dust contents Z dt iQ M and Z^ism- 111 
general, lowering the amount of dust raises the escape fraction, and the Lya line luminosity, 
and therefore increases the number of emitters above a fixed threshold, and vice versa. 
However, we find our results to be less sensitive to these changes than those in the star 
formation rate or in the clump covering factor, shown in Figure || In particular, when 
either Z^igm (dashed lines) or Z^ism (dotted lines) is varied between 0.1 and 1 Z , the 
corresponding range in the surface density of emitters is only about an order of magnitude 
at z ~ 4, and even less at higher and lower redshifts. On the other hand, the results 
are more sensitive to changes in the star formation efficiency (solid lines). For the range 
0.02 < < 0.2, the abundance varies by almost two orders of magnitude. Note that a 
low star formation efficiency results in an underestimate of the abundance that can not be 
compensated by changes in the other parameters (cf. Fig. ^). Finally, the long-dashed line 
in Fig. f| shows the effect of changing the cosmology to an open CDM model. As the figure 
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shows, around redshifts z = 3 — 5, the two cosmological models give almost identical results 
(although they differ at lower redshifts). We conclude that the parameters that are best 
constrained by the data are £*, -F cov , and e*. 

6. Conclusions 

We have studied the abundance of high-redshift Lya emitters by combining the 
Press-Schechter theory, describing the cosmological formation of halos, with a simple 
prescription for an inhomogeneous dust distribution and its dependence on the star 
formation process. We find that these simple models can explain reasonably well the 
recent observations and earlier non-detections of high-redshift Lya emitters. We find 
that it is possible to place significant constraints on the inhomogeneity of the spatial 
distribution of dust, parameterized by the clump covering factor -F cov , the star formation 
rate, parameterized by the duration of the Lya emitting phase and the overall star 
formation efficiency, e*. Our results show that in order to reproduce the surface density 
observed by Hu et al. (1998), F cov must be of order unity, t* must not exceed ~ 5 x 10 8 
yr, and the star formation efficiency must be at least ~ 5%. These numbers should be 
predicted by more complete and detailed models of galactic evolution, and will be useful 
discriminators between such models. 

Our models also predict that more Lya galaxies will be detectable, around the present 
flux threshold, upto redshifts as high as ~8, indicating that for a simple Press-Schechter 
cosmological model there can be significant Lya activity in the early phases of galaxy 
evolution. This is especially interesting in the light of recent discoveries of Lya galaxy 
candidates at redshifts z > 5. Three spectroscopically confirmed objects have been reported 
by Hu et al. (1998) at a redshift of 5.64, Weymann et al. (1998) at a redshift of 5.60, and 
Dey et al. (1998) at a redshift of 5.34. Clearly, the search for Lya emitters is beginning to 
enjoy considerable successes and these primeval objects may turn out to hold important 
clues both for galaxy evolution, as well as ISM physics. 

We thank Arif Babul for useful comments. ZH was supported by the DOE and the 
NASA grant NAG 5-7092 at Fermilab. M.S. has been supported by NASA through Hubble 
Fellowship grant # HF-01101.01-97A, awarded by the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc., for 
NASA under contract NAS 5-26555. 
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Fig. 1. — The escape fraction F csc of Lya photons from an inhomogeneous medium, as a 
function of dust content. The three curves correspond to three different values of the covering 
factor of opaque clumps, F cov =1, 5, and oo. 
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Fig. 2. — The derived evolution of the IGM dust content, Zd,iGM hi our models. The three 
curves correspond to three different values of i*. Since e* is kept fixed, the assumed star 
formation rate in each individual galaxy scales as 1/t*. 
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Fig. 3. — The surface density of Lja emitters in our standard model (solid lines) with fluxes 
above different values of the detection threshold. The two data points are taken from Hu et 
al. (1998). For the fixed threshold F = 1.5 x 1CT 17 erg cm" 1 s -1 , the dashed lines show how 
the surface density changes if the assumed star-formation rate is increased or decreased by 
a factor of 10. Similarly the dotted lines show the surface density when the covering factor 
is changed to F cov = 1, or oo. 
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Fig. 4. — Same as Fig. [| but now the solid lines show how the surface density changes with 
the star formation efficiency, e*; the dashed and dotted lines show the effect of changing the 
assumed IGM or ISM dust contents in the interval 0.1 — 1 Z ; and the long-dashed lines 
show the effect of changing the cosmology. 



